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Abstract 

Majoron models provide neutrino masses via the spontaneous 
breaking of a global U(l) symmetry. However, it may be argued 
that all global symmetries will be explicitly violated by gravita- 
tional effects. We show that it is possible to preserve most of 
the usual features of majoron models by invoking U(1)b-l to 
be a gauge symmetry and adding a second singlet scalar field. 
The majoron gets a small model dependent mass. The couplings 
of majorons to neutrinos may be of ordinary strength or may 
be made arbitrarily weak. We discuss the cosmological and as- 
trophysical consequences of majoron models in the context of a 
model dependent majoron mass and neutrino coupling. For an 
appropriate choice of parameters majorons can play the role of 
dark matter. 
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1 Introduction 



Majoron models Jl], @] were designed to accommodate a wide range of neutrino 
masses. In most theories, small neutrino masses are achieved via the see- 
saw mechanism ||. In grand unified models, seesaw neutrino masses are 
naturally of order 10~ 2 eV or less, but in majoron models the scale of lepton 
number violation is left as a free parameter, and so larger neutrino masses are 
not implausible. Another attractive feature of majoron models is that they 
contain a spontaneously broken global symmetry, complete with a Nambu- 
Goldstone boson, i.e., the majoron. This allows for neutrinos with mass 
greater than 30 eV to escape the constraint from the closure density of the 
Universe ||, since they may decay into a lighter neutrino and a majoron. 

Besides allowing for neutrino mass and neutrino decay, majoron models 
lead to a host of other cosmological, astrophysical, and laboratory conse- 
quences ||. Since the majoron couples to neutrinos there exist a number 
of constraints based on possible majoron production in neutrino rich envi- 
ronments; e.g. excess cooling of supernova cores by majoron emission || or, 
a majoron contribution to the energy density of the Universe which could 
destroy the successful predictions of big-bang nucleosynthesis 0. Couplings 
to other standard model particles arise through loops, which leads to con- 
straints based on the evolution of red giants || and horizontal branch stars. 
Majoron models also possess global U(l) cosmic strings, although we know 
of no important consequences of such relics. Finally, there are a number of 
laboratory experiments which are sensitive to majoron models. The most 
important examples are the width of the Z boson as measured at LEP M, 
which eliminates the simplest of triplet majoron models @; and the absence 



of 3-body neutrinoless double beta decay events |10[. We mention these 
features now because they illustrate the richness of majoron models, and 
they are all modified by our attempt to reconcile majoron models with short 
distance symmetry breaking induced by gravitational effects. 

A criticism of majoron models is that the global symmetry is not ne- 
cessitated by the known particle spectrum. Although there may be nothing 
wrong with this, it has been argued that quantum gravitational effects should 



explicitly violate all global symmetries [IT]. Assuming this to be true, a sim 



pie way to save the majoron models is to ensure that the global symmetry 
be accidental. This can be accomplished if we gauge U(1)b-l and choose 
the quantum numbers of the scalar fields to ensure that there are no gauge 
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invariant operators of dimension 4 or less which break the accidental symme- 
try. The higher dimensional operators which explicitly break the accidental 
U(l) will be suppressed by powers of M pi . 

A small explicit symmetry breaking has two immediate consequences. 
First, the majoron will acquire a small mass, which in turn has many astro- 
physical and cosmological implications. Second, in a cosmological setting, a 
U(l) symmetry that is broken both spontaneously and explicitly implies the 
formation of a network of strings and walls. The ability of this network to 
dissipate places constraints on the details of the model. 

In the usual majoron models there is a direct relationship between neu- 
trino mass and the coupling of neutrinos to majorons. In the present context 
there must be two scalar fields in order to break both the gauge symmetry 
and the accidental global symmetry. Both must have charge under B — L, but 
only one may couple to right-handed neutrinos to power the see-saw mech- 
anism. The majoron is predominantly the scalar with the smaller vacuum 
expectation value (VEV). If the smaller VEV is also associated with neutrino 
masses, then the linear relationship between neutrino mass and coupling to 
majorons is preserved; but if the larger VEV is associated with neutrino 
mass the coupling to majorons will be suppressed by a ratio of small to large 
VEVs. 

This paper gives a more in depth discussion of the issues just summarized. 
Section 2 reviews the singlet majoron model and its phenomenology. Sec- 
tion 3 discusses one way to build majoron models containing an automatic 
symmetry to protect against gravity. Section 4 discusses the cosmological 
constraints that arise from giving the majoron a mass. The possibility that 
majorons play the role of dark matter is also examined in this section. Sec- 
tion 5 is a discussion of the model constraints that arise from a consideration 
of the string-wall network. Sections 6 and 7 review the constraints from big 
bang nucleosynthesis and SN 1987A respectively. The last section contains 
a brief summary. 

Similar cosmological and astrophysical analyses have been performed for 
axion models with gravitational effects |T2|, |HJ , and for a broad class of mod- 
els where the mass of the Goldstone boson is due to instanton effects in an 
unspecified strongly interacting sector of the theory ||14|| . The situation for 
axion models differs from the majoron models presented here due to upper 
limits on the neutron electric dipole moment, which place severe constraints 
on the strength of explicit symmetry breaking due to gravitational effects. 
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These translate into severe constraints on the model, e.g. the lowest dimen- 
sional operator that is allowed from strong CP considerations is d = 10. No 
similar constraint exists for majoron models, and so d = 5 is viable. 

The phenomenology of textures [15| and neutrino masses [16| have also 
been reexamined in the context of planck suppressed symmetry breaking. 
While this work was in preparation we became aware of a paper by Akhmedov 
et al. [17] which considers some of the same issues as the present paper for 



majoron physics. 



2 The Minimal Majoron Model 

The simplest majoron model U is built by adding SU(2)l x U(1)y singlet 
neutrinos (N') Q as well as a singlet scalar (</>) to the spectrum of the standard 
model. The singlet scalar couples only to the singlet neutrinos because it 
carries lepton number L = —2. The relevant part of the Lagrangian which 
carries these extra fields is given by 

L new = ^N!i?Nl + V(H,<i>) 

+ ^\WNl(l + l5 )N;<j ) +^ ) ui(l + l5 )N;H + h.c^ , (1) 

where H denotes the standard model Higgs doublet and i,j are generation 
indices. We use four component majorona spinors (N c = N, v c = v), and 
assume, for simplicity, that the Yukawa coupling matrices \^ contain 
no CP violating phases. The Higgs potential can be written as 

V = k x [HH^ - v 2 H } 2 + k 2 [HH ] - v 2 H }[4>4>* - v 2 } + k 3 [00* - v 2 } 2 . (2) 

If 4kiK3 > then (H) = vh and (0) = v. Utilizing the non-linear realiza- 
tion of the global symmetry, the scalar may be written as 

<f) = (p + v)exp(ix/v), (3) 

where the majoron will be \i an d the real scalar piece p will be massive. We 
can eliminate the majoron field from the mass matrix by transforming the 

*We use primes to denote the interaction basis for the neutrinos and leave off the primes 
to denote the mass basis. 
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neutrino fields; N- — > exp(—ixj5/2v)Nl and v[ — > exp(ix7 5 /2u)z/ 4 / ; leaving 



+ dpx&fl^i ~ Nh^M/(Av), (4) 



where 



wiu.ti = A§ M ij = A§ 5 «, (5) 

the subscript .D indicates a Dirac mass term coupling z/ and N', and we have 
not written the coupling of neutrinos to p particles since we are not interested 
in theories where the p is light. This will yield the following see-saw mass 
matrix 

L mass = 1/2 ( v> N')m( V n , J (6) 

where 

*-( A %*)■ (7) 

Upon diagonalization this mass matrix gives the usual see-saw spectrum of 
three light neutrinos (u) with masses m v ~ (m 2 D /M) and three heavy neu- 
trinos (N) with masses 0{M). 

It was pointed out long ago that the mass matrix and the majoron 



coupling matrix for the light neutrinos will commute up to 0((i71d/M) 4 )). 
Therefore, unless we make the heavy scale v < 0(100 GeV) the lifetime 
of any neutrino v with mass greater than 30 eV will be too small to avoid 
the constraint from the age of the universe, not to mention the more severe 
constraint on the lifetime from structure formation ||19|| . As a result, it has 
been generally accepted that additional Higgs fields must be added in order 
to get viable models for neutrino decay P0 |. However, it has been recently 



pointed out [^TJ that the result in reference Jl8[ does not hold in general. 
Radiative corrections will generate finite pieces of 0(m 2 D /(8ir 2 M 2 )) to the 
off diagonal part of the light neutrino couplings which will not be the same 
as the corrections to the off-diagonal terms of the mass matrixF]. The light 

2 It is also possible to get rapid decay if one chooses the global symmetry to be a linear 
combination of the lepton numbers (L e ^ ;T ) as opposed to just L p2j . Rapid decay into 
three neutrinos is possible in models with Dirac neutrino masses and extra singlet scalars 
which do not acquire VEVs [f23f. 



5 



neutrino Lagrangian may then be written with derivative couplings as 

1 



2 



+ \ Vil^l^i + J2 P rf{9v,i3 + 9a,ijlh)Vj ] , (8) 



I.) 



where g v and g a are of order 0(m 2 /(8ir 2 M 2 )) and there may be vector off- 
diagonal couplings, i.e. g V) a = 0. The decay rate for the heaviest of the light 
neutrinos into a particular channel is then given by 



Tfly I Tfly \ 2 2 2 



64lT V V 



(9vAj + g a ,ij)> ( 9 ) 



where we have neglected the mass of the final state neutrino. 

In the singlet majoron model, coupling to other fermions arises through an 
intermediate Z boson and a neutrino loop, see Fig. |l[ The electron coupling 
in pseudoscalar form is then 

G F m e m v 

W ^6^- (10) 

The most important constraint on this coupling arises from consideration of 
the cooling of red giant stars 0, g eex > 1CT 13 . 

3 The Singlet Majoron Model with Explicit 
Symmetry Breaking 

Assuming that short distance physics only respects local gauge symmetries, 
the only way to have Nambu-Goldstone bosons appear in the spectrum is to 
ensure that there is an accidental (clearly a misnomer in this case) symmetry. 
There is no way to prevent operators which break the accidental symmetry 
from appearing at some level, but we can arrange for these operators to have 
dimension greater than four. In this case the symmetry breaking operators 
will be suppressed by powers of M pi and a small mass is induced for the 
Nambu-Goldstone boson. 

To accomplish this for majoron models, we gauge U{1)b-l and include 
in the spectrum two singlet scalar fields. The first, <fii, plays a similar role to 
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(j) in the singlet model of section 2. This field has B — L quantum number 
qi = 2, allowing it to have a Yukawa coupling with the singlet neutrinos. 
The U{V)b-l charge of the additional scalar field, 2 , is chosen such that 
the lowest dimensional B — L invariant operator which breaks the accidental 
U(l) has dimension greater than four. To leave matters as general as possible 
we allow <p2 to have fractional B — L charge. As such, there will be several 
different models possible for a given dimension; for example, g 2 = 1/2, 4/3, 
3, or 8, will result in a dimension 5 operator suppressed by one power of M v \. 
As examples of g 2 assignments that do not work, g 2 = 1 or 4 allow a planck 
scale cubic term in the lagrangian and the global symmetry would disappear 
entirely from the low energy lagrangian. 

We write the lowest dimensional symmetry breaking term as 

Ls.b = g 9rav <t>r<t>?*/Mir A + h.c, (11) 

where g grav is the dimensionless strength of the operator, d = n\ + n 2 is the 
dimension of the symmetry breaking operator, and ni and n 2 are integers. 
Gauge invariance imposes the constraint 

ni<?i - n 2 q 2 = 0, (12) 

and to ensure that d > 4 we require 

aqi — bq 2 ^ : if a + b < 4; a, b G integers. (13) 

By choosing an appropriate Higgs potential both singlets will receive vac- 
uum expectation values, v\ and v 2 respectively. Given that B — L is a local 
gauge symmetry there will be a lower bound on the Z' mass of 0(300 GeV) 
which, for reasonable U(1)b-l gauge couplings, will also be a lower bound 
on the larger of the two singlet VEVs. Had we chosen to gauge an anomaly 
free combination of lepton numbers this constraint would be relaxed . 

Upon symmetry breaking, one combination of the phases of the singlets 
becomes the longitudinal part of the Z' gauge boson, while the orthogonal 
combination becomes the majoron, 

The diagonal majoron couplings to the light neutrinos are given in the pseu- 
dovector form as 
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where the first factor accounts for the mixing of 0i and 02 in Eq. and for 
convenience, we define 

vl = f v * (16) 
e// qlvl + q 2 2 v 2 2 

The explicit symmetry breaking will induce a mass for the majoron, 

m no 

2_ 9 VV n 2 n 1 

The coupling of other fermions to the majoron arise from the same dia- 
gram as in the simple model of section 2. At first it might appear that one 
could replace the Z boson by a Z' boson, while at the same time replacing 
the light neutrino loop by a loop with a heavy singlet neutrino. If such dia- 
grams existed they would dominate over the Z exchange by a ratio of heavy 
to light neutrino masses, however, one can show that such contributions van- 
ish. First, the \ — Z' mixing (through a heavy singlet neutrino loop) cannot 
lead to /75/x vertex since B — L is a vectorial symmetry. There will be 
Z — Z' mixing at one-loop via the exchange of charged fermions and, since 
Z coupling to matter is parity violating, /75/x may be induced as a two- 
loop effect. However, such diagrams also give vanishing contribution. This 
is best seen by working in Landau gauge, where the gauge boson propagator 
is purely transverse. There is no x ~ or X ~ Z mixing in this gauge, so 
only the unphysical Higgs may contribute; but this scalar is necessarily an 
S'?7(2)i-singlet (since it has to couple to the heavy singlet neutrino) and so 
has no coupling to the standard model fermions. 

The constraint from coupling majorons to ordinary matter is therefore 
the same as in the simplest model, with the exception that even this mild 
constraint does not hold if m x > 10 keV, since the majorons will be too heavy 
to be produced in red-giant cores. 



4 Majoron contributions to the energy den- 
sity 

In theories where a global symmetry is broken spontaneously as well as ex- 
plicitly by a term with a relatively small coupling, there are two contributions 
to the energy density in Goldstone bosons, the thermal production due to 
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interactions with the other particles in the medium and a coherent oscillation 
of the Goldstone field if the initial angle of that field does not align with the 
direction of explicit symmetry breaking. The coherent contribution may be 
enhanced by the presence of strings and walls associated with spontaneous 
symmetry breaking, or may be suppressed by an inflationary epoch. Both 
terms may be alleviated if the Goldstone bosons decay rapidly. 

We begin by summarizing the cosmological history of the model. For this 
task it is convenient to refer to the scalar fields by the magnitude of their 
VEVs, in which case we will refer to the smaller VEV (and its field) by the 
subscript a and the larger VEV by the subscript &. It may be that a = 1, 
indicating the scalar field that couples to neutrinos, or that a = 2, depending 
upon the scenario. Then, the larger of the VEVs, u&, is established first, the 
associated scalar is "eaten" by the U{1)b-l gauge field, and the Z' gets a 
mass. At this time a set of local U(l) gauge strings will form. 

When the remaining global U(l) symmetry is broken, v a is established 
and the majoron becomes a light degree of freedom. Breaking the global 
symmetry results in a set of global cosmic strings. Additionally, the local 
gauge strings present from the breaking of U(1)b-l ma y become catalysis 
sites for another set of global strings. 

The majoron has a small mass (Eq. [17D which is present immediately at 
the breaking of <p a . This mass becomes dynamically important at a time 
labeled by the temperature T m , defined by 

m x = 3H(T m ), (18) 

where H is the expansion rate of the Universe. It is possible that m x > 
H at the time when v a is established, in which case T m rs v a . In either 
case, two things occur at T m . First, one may expect a contribution to the 
energy density from a non-thermal population of long-wavelength coherent 
majorons. Second, the global cosmic strings become the edges of dynamically 
important domain walls, whose thickness is of order l/m x . Depending upon 
details of the model these string-wall networks may or may not go away. 

At first the majoron field should be thermally coupled to the plasma 
through scattering processes, but at some time (labeled by the temperature 
T d ) the majorons will decouple. It is possible that decays and inverse decays 
may recouple the majorons to the plasma at a later date. The rate of inverse 
decays in a relativistic plasma is ~ T x m x /T, whereas the expansion rate 
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IS ~ T 2 /M pl . Evidentially, inverse decays can bring a species back into 
equilibrium at low temperatures - depending on the interaction strength. 
If the majoron-neutrino coupling is weak then full equilibrium will not be 
achieved until the majorons are non-relativistic - at which time their energy 
density will be dumped into neutrinos but back reactions will be unimportant. 
Of course, this may occur after the present epoch. 



4.1 Thermal contribution to the energy density 

If the majorons are lighter than a few eV their thermal contribution to the 
energy density will not be significant, but a larger mass is certainly possible 
in the present context. The energy density in a thermal relic population of 
majorons is 

Pth = n th m x = ^-T^m x , (19) 



where T x is the 'temperature' of the majorons. Generally, T x will differ 
from the photon temperature by a factor R x , in which case the thermal relic 
density can be given in terms of the critical density for closure, p cr u, as 

Peru keV /ig 

where ho is the expansion rate in units of 100 km sec -1 Mpc^ 1 , and we used a 
microwave background temperature of 2.74 °K. If the Universe has expanded 
adiabatically since the majorons decoupled, and light neutrinos are the only 
non-photonic species present today then 




i$=pO = 91 ' 7 + ' £ + ^ U5.8xlO- 2 , (21) 

' ; V//- • .'/- //• / 



where g 1 — 2, g e — 7/2, g u = 21/4, and g* is the total number of relativis- 
ts degrees of freedom other than majorons at the time that the majorons 
decouple. 

As will be shown below, for most scenarios decoupling may be expected to 
occur soon after a develops a VEV. At first the self couplings of the majorons 
and p particles and their coupling to singlet neutrinos may be expected to 
keep the majorons in equilibrium. If v a and vj, are not too dissimilar the 
B — L gauge interactions will couple the singlet neutrinos, majorons and p 
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particles to the ordinary degrees of freedom in the plasma. However, all these 
particles, except for the majoron, presumably have masses of order v a or vt,, 
and will soon freeze out of the primordial plasma. After they freeze out the 
only way to avoid decoupling is through the light neutrinos. The most direct 
reaction to maintain equilibrium is vv <-> XX through t-channel exchange of 
a neutrino. The cross section for this process is given by 

ml fq 2 v e ff 



This result applies in the limit where the p particles are much heavier than the 
momentum exchange in the reaction. The last factor arises from projecting 
the majoron onto the scalar field that couples to the neutrinos. Allowing this 
factor to be of order unity, the decoupling point is given by 



Unless f i is quite small and the neutrino masses are large the light neutrinos 
cannot keep the majorons in equilibrium after the phase transition. The 
converse of this is that the constraint on majoron models from big bang 
nucleosynthesis is fairly weak (see Section 6). Thus, except for exceptional 
circumstances, it is reasonable to assume that the majorons decouple soon 
after the phase transition. 

Eq. |20] is valid only so long as the majorons are stable. However, if 
majorons have a mass then they may decay with a rate 

r* = ^I>i. (24) 

where 

= (25) 

Vl Vi 

is the diagonal coupling to neutrino i in its pseudoscalar form. The sum is 
over neutrino species lighter than the majoron. Off diagonal couplings are 
presumably less important. The decay rate depends on the light neutrino 
masses, which arise through the seesaw mechanism. Phenomenological con- 
straints on neutrino masses will impose constraints on the majoron decay 
rate. 
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Some general remarks, a) If neutrino masses are small enough, then 
majorons are stable for all practical purposes, b) If the majorons decay well 
after they become non-relativistic, the decay neutrinos will contribute to the 
energy density in radiation, and this can produce a conflict between theories 
of galaxy formation and limits on the fluctuations in the temperature of 
the cosmic background radiation. In the extreme case the density of decay 
products can overdose the Universe, c) If the majorons decay too quickly 
then inverse decays may recouple the majoron field to the plasma before 
nucleosynthesis begins at a time of about 1 sec, or a temperature Tbbn ~ 
1 MeV. In this case an excess of 4 He will be produced (see section 6). d) 
If majorons are more massive than about 10 MeV then their decay before 
nucleosynthesis (more specifically - before the neutrinos decouple) poses no 
significant cosmological problems other than a potential increase in entropy 
and consequent dilution of baryon number. 

Let us investigate the possibility that relic majorons form warm dark mat- 
ter. As an example, consider the case where the two VEVs are comparable, 
V\ — v-i — v. We vary q 2 , the B — L charge of 02, so as to get the appropriate 
dimension operator. To get the majoron to be warm dark matter we want 
the mass of the majoron to be of order a keV. Using Eq. |T7| and assuming a 
dimensionless strength of order unity for the explicit symmetry breaking, we 
find a relation between v and the dimension of the operator, 

log 10 (Vl GeV) ~ 19( ^~_ 4) 2 ~ 12 . (26) 

The decay rate will be given by Eq. If we scale to a neutrino mass 
of 10 eV, keep the majoron mass at a keV, and drop the numerical factor 
(q2V e ff/vi) 2 then the majoron lifetime is 

o 1A J f V/iOeV^ 2 /keV\ 

r *~ 3xl ° y (— ) U:) «*■ < 27) 

If d = 5 then v ~ 200 GeV and the lifetime is around 10 4 seconds. Majorons 
have mass but decay quickly so they do not constitute dark matter. For 
keV majorons, their decay rate lies in an acceptable window which poses 
no problems for galaxy formation or for nucleosynthesis. If d — 6 then 
v ~ 3 x 10 6 GeV and the lifetime is around 3 x 10 12 sec. This decay rate 
is troublesome for galaxy formation, but a modest decrease in the heaviest 
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neutrino mass would allow the majoron to be sufficiently stable to be dark 
matter. For d = 7 we find v = 10 9 GeV and lifetimes longer than the 
age of the Universe, so majorons can play the role of warm dark matter 
in such models. Higher d also admit a value of Q x ~ 1 in stable thermal 
relics. Clearly, we could vary the ratio of the two VEVs, the dimensionless 
strength of the gravitational operator, etc., providing even more flexibility 
to the models. 

Let us look at the d = 5 little more closely. This can be achieved 

by qi = 2 and q 2 = 1/2, 4/3, 3, or 8. Let us concentrate on q 2 = 1/2, 
and still take the two VEVs to be the same. Assuming that both VEVs 
develop before the standard model breaks from SU{2) x U(1)y down to 
U(l) em , so that g* = 67|, we find R x = 0.058. For a Hubble parameter of 
h = 0.5 (where h is H in units of 100 km/ sec/ Mpc) and a majoron mass of 
53gll 2 av keV(u/l TeV) 3 / 2 this gives fi x = 242#V 2 >/1 TeV) 3 / 2 , if the 
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majorons were stable. However, with these parameters and a 30 eV neutrino 
mass, the lifetime of the majorons is of order 10 4 seconds. 

Coincident ally, the temperature of the Universe at 10 4 seconds is of order 
10 keV, so the majorons would be mildly non-relativistic when they decay. 
The decay neutrinos would be slightly more energetic than the blackbody 
neutrinos, but not so much as to pose a problem for galaxy formation. How- 
ever, as Akhmedov, et al. [|17| point out, if one assumes an upper limit on the 



neutrino mass of roughly 30 eV, then this argument yields an upper limit to 
the VEV, v ^ 10 TeV. Both m x and r x increase with v and so the present 
density in relativistic decay neutrinos also increases; which eventually creates 
a problem for galaxy formation. In that analysis there is only one VEV, but 
in our models there is more flexibility. By separating v% and and choosing 
different B — L charges one can maintain a short majoron lifetime while in- 
creasing m x . For example, in the toy case we chose q% — 2 and q 2 = 1/2 so 
that varying v\ and v 2 independently (but keeping v\ < v 2 , so that v e ff ~ Vi) 
gives 

- 1 / 2 2 

m x ss 1 1/2 2 28 



8tt / vt \ 2 



m x \m L 



(29) 

As v 2 is raised the majoron mass increases but the lifetime shortens since 
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the neutrino couplings haven't been altered. The energy in relativistic de- 
cay products increases more slowly. There is still a prohibited region from 
BBN considerations, 0.1 MeV 5; m x £ 10 MeV, but generally more massive 
majorons are still viable in a model with two VEVs. 

Another interesting case occurs when the heaviest neutrino has a mass 



chosen to solve the solar neutrino problem via the MSW mechanism [£5 
m v ~ 10~ 2 ' 5 eV. The majoron lifetime would be around 10 12 seconds in 
such a model. This value is somewhat uncomfortable - the decay products 
would be relativistic and the Universe would be dominated by radiation for a 
period late in its history, making a consistent theory of galaxy formation more 
difficult. To rectify this, it would be desirable to make the majoron lighter 
and longer lived than the age of the Universe. This could be accomplished by 
increasing v% while holding t> 2 fixed and the neutrino masses constant. This 
would decrease the majoron mass and their coupling to neutrinos. 

Finally, one of the main motivations for majoron models is to allow for 
neutrino masses in excess of 30 eV. The presence of the massless majoron 
allows the heavy neutrino to decay into a lighter neutrino and a majoron with 
a short lifetime. The issue of off diagonal couplings may be addressed by 
invoking radiative loops, or enlarging the model to include extra scalar fields 
so that the mass and coupling matrices are not nearly proportional. With 
the consideration of gravitational effects one may expect the majoron to get 
a mass, perhaps in the keV region. Neutrinos heavier than the majoron may 
still decay into a majoron and a lighter neutrino, but neutrinos less massive 
than the majoron, and yet more massive than 30 eV have no obvious decay 
channel and would present a cosmological problem. 

4.2 Coherent Oscillations 

The small explicit breaking of a spontaneously broken global symmetry leads 
to a coherent oscillation of the Goldstone field around its minimum. This 



phenomenon is well known in axion models |26| and leads to a constraint 
on the axion decay parameter, f a . If gravitational effects explicitly break 
the Peccei-Quinn symmetry in an axion model, then the energy in coherent 
oscillations is smaller than that in thermal axions unless the explicit symme- 
try breaking operator has d > 10 [[T^, |I3" |. The models considered here are 



similar except for the presence two VEVs, but we will show that the second 
VEV can only make the coherent oscillations less important. 
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The majoron field is static until the mass becomes greater than the ex- 
pansion rate of the Universe, m x = 3H(T m ). The energy density in majorons 
at that time is 

Pmis = 7; m lx 2 92qiVeff ; (30) 
2 A n 2 ni 

where xo is th e initial misalignment angle, presumably of order unity. These 
coherent oscillations have a wavelength which is greater than their Compton 
wavelength and the energy density scales as non-relativistic matter. Taking 
into account the expansion of the universe from T m until today yields 

Q^TxlO-^^^ML, (31) 
x h 2 GeV 5/2 n 2 nx v ; 

where g m is the number of equilibrated degrees of freedom at T m . 

To compare the thermal relic and coherent energy densities we compare 
their contributions to the majoron number density at the temperature T m . 
Dropping all dimensionless factors, the ratio is given by 

n m V 2 ( \ ( 10 - d )/ 4 /„ x n b /4 

n mis m X U eff V a \ ( v a \ 



Rn x = — * 3/2 3/2 ^ P ( 32 ) 

nth m 3 x /2 Mp{ 2 \MpiJ \v b J 

where we have used T m w (m x M p ;) 1//2 in the first relation; and Eq. |l^for the 
majoron mass and Eq. |16] for v e ff in the second. Note that the result depends 
only on the relative magnitude of the VEVs and not on which one couples 
to neutrinos. The power (10 — d)/4 is the same as in the axion result. The 
extra ratio of VEVs reflects the fact that a large Vb will increase m x without 
affecting f e //, but this can only reduce R n . Finally, for d — 5 (and maybe 
for d — 6) T m w v e ff, i.e. the majoron mass is greater than the expansion 
rate at the moment of the phase transition. In this case R nx ~ m x /v e ff, and, 
again, misalignment is unimportant. Thus, the misalignment energy is only 
important if the dimension of the gravitational operator is d > 10. 

Eq. [31] may need to be modified to account for the majorons that arise 
from cosmic strings. Following the arguments made by Davis and Shel- 



lard |23] for axion models, the presence of global strings may be expected to 
increase the number of majorons in the coherent field by perhaps an order of 
magnitude due, primarily, to the logarithmic divergence in the string tension 
of global strings. However, this conclusion depends on the spectrum of Gold- 



stone bosons radiated from cosmic string loops at the time when T = T m 2c 
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A conclusion to the argument awaits a definitive numerical analysis of the 
radiation spectrum. 

The derivation of Eq. |3T] also assumed that the initial misalignment angle 
was of order one. If there were an inflationary epoch between the breaking of 
the global U(l) symmetry and T m then the majoron field could be correlated 
over horizon size regions and could, by chance, be near its minimum. This 
would yield a reduction in the misalignment energy density from the estimate 
in Eq. [30] due to a suppression of xo- The reduction is limited by quantum 
fluctuations in the majoron field during the inflationary epoch to a factor 
of order (Hi/v e ff) 2 , where Hj is the expansion rate during inflation. One 
must also check that isocurvature density perturbations induced by these 
fluctuations are not in conflict with observations of the cosmic microwave 
background radiation. 

One may wonder if the presence of a small explicit symmetry breaking 
can cause the initial angle of the majoron field to lie near the minimum 
of the majoron potential, i.e. xo ~ 0. The same issue is relevant for the 
question of the formation of majoron strings and domain walls (see below). 
We believe the answer is no, the angular minimum is not strongly preferred. 
The essential point is that the radial gradient of the potential will be much 
greater than the angular gradient caused by the explicit symmetry breaking, 
see Fig. 0. Once a VEV has been established, it will roll radially outward 
much faster than it can slip sideways. The one possibility seems to be that the 
tilt causes the VEV to be established in a direction aligned with the explicit 
symmetry breaking, so that the radial roll of <fi takes one directly to the global 
minimum. To study this, we consider a first order phase transition which 
proceeds through the formation of bubbles, which may form by quantum 
tunneling or by thermal fluctuations. We discuss the tunneling case first and 
comment on thermal fluctuations later. 

The rate of bubble formation depends upon the action to form a bubble 
of zero energy, which in turn depends upon the amount of supercooling and 
the shape of the potential. Dimensional analysis suggests that the action for 
a critical size bubble is 

A c ~ R*AV ~ R 3 c a ~ (33) 

where R c ~ a/AV is the radius of a critical size bubble, AV^ is the potential 
difference between the false and true vacua, a ~ v 2 /l ~ AVI ~ v(AV) 1 ^ 2 is 
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the surface tension on the bubble, I ~ vj (AV) 1 / 2 is the surface thickness, and 
v is the magnitude of the VEV to be established after the phase transition. 
The formation rate per unit volume scales as 

T c ~ R-'e-^. (34) 

This rate should be compared to the expansion rate of a horizon volume, 
~ H A . Now, even though the prefactors may have ratios of large numbers 
such as Mpi/T this comparison will be dominated by the exponential. When 
supercooling first begins AV will be quite small, the action for bubble for- 
mation will be large, and bubble formation will be slow compared to the 
expansion rate. However, as the Universe continues to expand and super- 
cool, AV will increase and A c will decrease. Eventually, A c will be of order a 
few hundred and bubble formation will become rapid. The significant point 
for us is that by the time this happens the tilt of the potential will be rather 
small, e <C AV. As a result the action to form a critical bubble with the 
majoron field at a local maximum will be of order R^(e + AV) whereas the 
action to form a bubble at the global minimum of the potential energy is 
just i?^AV. Since bubble formation proceeds when A c ~ a few hundred, we 
require that tj AV ~ 10~ 3 in order to significantly bias bubble formation. 
Since AV is of order Xv 4 , this requires a very small self coupling for the scalar 
fields, which in the context of a renormalizable U(l) local gauge theory is not 
natural. We conclude that for the case of quantum tunneling the majoron 
field is randomly aligned with respect to the explicit symmetry breaking. 

The thermal fluctuation case is similar. Instead of considering exp(— A c ), 
one considers exp(— /3F C ), where F c is the free energy of a thermal fluctuation 
forming a bubble for which the energy is maximal, dE/dR = 0. Which 
process dominates depends on the ratio f3F c /A C) which in turn depends on the 
details of the theory. Generally, the transition for a weakly self-coupled scalar 
field should be dominated by thermal fluctuations, whereas for a strongly 
coupled field the tunneling process should win out. In either case, the explicit 
breaking due to gravity should not be important. Therefore, for a first order 
phase transition, even though the majoron mass may be greater than the 
expansion rate at the time of symmetry breaking, there will not be any 
suppression of the coherent energy density due to a physical alignment of the 
majoron field along a preferred direction, and Eq. |31] may be used. 

As with the thermal majorons there is the possibility of majoron decay; 
however, if we confine ourselves to the case where the coherent energy is more 
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important than the thermal energy we can show that majoron decay is not a 
relevant issue. By using the ratio of coherent to thermal majoron densities, 
Rn , given in Eq. we can express the majoron mass as 



m v = 




£ GeV. 



(35) 



Using this form in Eq. [H] gives 



V eff 



10 8 WX /4 GeV. 



(36) 



To maximize the decay rate we take the neutrino mass to be comparable to 
the majoron mass, take v\ = v e ff, and find 



Decreasing the lifetime requires either a large value of Q x with a correspond- 
ing unacceptable density of decay neutrinos, or requires making a thermal 
density of majorons in excess of the coherent density. 

5 Majoron Strings and Walls 

In the models presented here there are several types of topological defects 
which may occur. When U(1)b-l breaks, a set of local cosmic strings forms 
with string tension fn, ~ v\. We call these type-6 strings since they have a 
defect in the field 0&. Later, when T ~ v a , the second U(l) breaks and type-a 
global strings form with \x a ~ v\. At the same time, type-6 strings develop 
a defect in v a with winding number adjusted to minimize the energy in the 
angular gradients of the fields. Since there is explicit symmetry breaking in 
our models, at the time T m one expects domain walls to form with surface 
tension a ~ m x v 2 e ^. These domain walls may exist as isolated surfaces or be 
bounded by strings. If m x > H(v a ) these walls will exist from the moment 
that majorons exist, but otherwise they form at a later time. 

The fate of the string-wall network depends upon how many walls are 
connected to the various types of strings. If either the type-a or the type-fe 
strings come with a single wall attached, then the system can dissipate by self 



r. 



x 



> 




sec. 



(37) 
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intersections |29j. If, however, both types of strings have a number of walls 
N w > 1, then the strings cannot annihilate and the walls may eventually 
dominate the energy density of the universe. The energy density in a single 
wall stretching across today's horizon is 

~ WW vl ~ 10 "GeV^ (38) 

A single wall stretching across the horizon would be likely to contribute to the 
quadrupole moment of the microwave background, and so it seems reasonable 
to require Vt wa u £ 1CT 5 . The resulting constraint, m x v 2 e ^ < 1CT 9 GeV 3 is much 
stronger than that which follows from the coherent energy density (Eq. |3l|) , 
and somewhat stronger than the constraint from thermal relics, depending on 
the magnitude of f e //- For example, it is violated by ~ 7 orders of magnitude 
for the toy model used in presentation of the thermal relics. If majorons are 
to make a significant contribution to Q, the model must be defined so that 
N w = 1 for some set of strings. On the other hand, it is not difficult to 
suppress m x by sufficient powers of M p i in order that majoron walls pose no 
problem. 

Constructing models where N w = 1 for either type-a or type-6 strings 
may be accomplished by an appropriate choice of the B — L charges q a and 
q b . The problem is identical to that considered in reference fl3| for a class 



of axion models, so we summarize that analysis. The energy density in the 
majoron field around a string is proportional to 




~ {qbW a + q a w b ) ~ {n b w b + n a w a ) , (39) 

where w a and w b are the winding numbers for <ft a and <p b respectively. Recall 
that n a /n b = qb/q a , with both n a and chosen to be integers with no 
common factors other than 1. At the same time, the number of walls is given 
by the number of maxima in the gravitationally induced potential (Eq. [TT| ), 

N w = n b w b + n a w a . (40) 

Therefore, minimizing the field energy reduces the number of walls to a 
minimum as well. 

For minimal type-6 strings we have w b = 1, and w a is chosen to minimize 
the gradient energy; while for type-a strings w b = and w a = 1. There are 
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three possibilities for achieving strings with but one wall, a) If n& = 1, type-6 
strings have w a = and one wall, b) Similarly, if n a — 1, then type-a strings 
have a single wall, c) It is possible for n a ^ 1 and n b ^ 1 but to still have 
that w a which minimizes Eq. |39] result in a single wall; e.g. if q a = 2 and 
qb = 6/5, then n b = 5, n a = 3, and w a = 2. 

Even if the network can dissipate, there is an additional constraint that 
this happens quickly enough. The system cannot dissipate until the energy 
in walls is comparable to the energy in strings. If the string energy dom- 
inates, then small loops of string that intersect walls result in holes which 
are subcritical in size; i.e. they do not expand to eat up the wall, but rather 
the hole shrinks and fills in. In order to get dissipation the critical size hole 
must be less than the horizon size. For type-6 strings, this occurs when 
(vb/v a ) 2 (H/m x ) > 1. This may be later than the time when the majoron 
mass first becomes effective due to the difference in the magnitude of the 
VEVs that control the characteristics of the walls and the strings. If this 
time occurs after the time of matter-radiation decoupling, and v b ~ 10 16 GeV, 
then the strings still present at that time will cause excess fluctuations in the 
microwave background. 

There are additional constraints that arise from the decay products of 
the strings and walls. For example, there are the majorons created in the 
dissipation of the type-a strings, which strengthen the constraint from the 
misalignment energy, Eq. [IT} Other examples concern the type-6 strings. 
After they first form and before they get dressed up with a majoronic co- 
coon at T ~ v a , a scaling solution is established and presumably maintained 
by the radiation of gravitons. Requiring that the graviton energy not con- 
tribute too much density during nucleosynthesis provides a constraint on Vb, 
Vbln(vb/v a ) ~ 10 18 GeV [[KJ. Similarly, if n a = 1, type-6 strings develop no 
walls at all for w a = —rib. In this case the number density of type-6 strings is 
always determined by a scaling solution. Requiring that density fluctuations 
induced by these strings do not cause excess fluctuations in the microwave 
background implies i>& ^ 10 16 GeV. 

6 Thermal Majorons and Nucleosynthesis 

If majorons are in thermal equilibrium at T = 1 MeV they will contribute 
to the energy density an amount equivalent to 4/7 of a neutrino species. 
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However, Walker et al. [31] argue that based on a determination of the pri- 
mordial Helium abundance the excess energy density cannot exceed the stan- 
dard cosmology by more than 0.3 species, i.e. AN BBN < 0.3. As pointed 
out earlier there are two ways for neutrinos to come into equilibrium - either 
through scattering events or through decays and inverse decays. Bertolini 
and Steigman j7| have examined the scattering case for the simple singlet 
model discussed in section 2. Using the cross-section of Eq. |22| we adapt 
their result to the model in section 3 and find 

^vf— )7^V< 4.6x10- (41) 
MeV V v 1 J V vx J v 1 



where the last factor accounts for the fact that the majoron does not neces- 
sarily couple to the neutrinos. If Eq. fy] is satisfied the majorons freezeout 
before T ~ 100 MeV and, as a result, the entropy released when the pions 
and muons freeze out dilutes the majorons and their contribution to the den- 
sity is less than 0.3 neutrino species. If the majorons freeze out later they 
share in the entropy release and the Helium abundance does not come out 
right. 

The other constraint comes from decays and inverse decays p2| . As a 
safe constraint, we require that inverse decays come into equilibrium after 
the temperature drops below 5 MeV. If the majorons come into equilibrium 
before that, then we can be sure that they would share in the entropy of the 
Universe before nucleosynthesis begins and before the neutrinos themselves 
freezeout. If the inverse decays do not become effective until later they will 
only share the neutrino energy density. The effects on Helium production 
are more difficult to calculate (see, for example, Dodelson and Turner []33 



for a discussion of the subtleties in phase space distributions during neutrino 
freezeout). We therefore take as an approximate condition 

r -37Tkv s H{b MeV >' < 42) 

where the subscript , indicates either neutrino or majoron depending upon 
which is heavier, Tj is the decay rate at rest, and the second factor accounts 
for the Lorentz dilation of the decay of relativistic particles. The factor of 
3 multiplying the temperature crudely accounts for the average energy of a 
massless species. We ignore cases where two particles are nearly degenerate 
in mass, and we take no account of quantum statistics (Pauli blocking, Bose 
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enhancement) in determining our criteria. Evaluating Eq. [|2]at Td ec = 5 MeV 
when the number of degrees of freedom is g e ff = 43/4, 



MeV sec . , . . 

> 1.3 x KT 3 . (43) 



mi Ti 



We first consider neutrino decay in the simple singlet model of section 2. 
Assuming the neutrino mass is less than 5 MeV, the constraint is 



where 9 = Jga,ij + 9v,ij ( see Eq. 0) is an off-diagonal coupling strength ex- 
pected to be much less than one. If 9 and q%v e jjjv\ were equal to 1, the 
decay constraint is stronger than the scattering constraint for v > 14 GeV, 
but for 9 = 10 -3 they are comparable at v — 140 GeV. For m u < 30 eV and 
v > 100 GeV neither constraint is restrictive, nor will either constraint be- 
come more strict if we dilute the neutrino- majoron coupling by adding extra 
scalars. 

With explicit symmetry breaking the majoron gets a mass and we must 
consider the reactions, \ <-> vv. These reactions are more interesting than 
neutrino decay since a) they can proceed through the diagonal neutrino ma- 
joron coupling and b) the majoron mass can be larger than the neutrino mass 
and this allows for a faster decay. As an example, consider our toy model; 
with v\ = t>2 = v, qi = 2, q 2 = 1/2, and g gra v = 1, the constraint is 

ssb) <sv 5 5L (45) 

If we confine ourselves to m v < 30 eV then the constraint on v is quite weak, 
v < 10 10 GeV. This is a fictitious constraint - the majoron mass increases 
above 10 MeV for v»5x 10 4 , and for masses above this our analysis does 
not apply. In Eq. ^S|the dependence on v arises due to the t> 3//2 dependence of 
m x . In other models this dependence will differ, and so will the constraints, 
but generally for neutrino masses below 30 eV big bang nucleosynthesis does 
not pose any serious constraints. 
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7 Supernova Constraints 



Since the SN 1987A neutrino pulse had its expected duration (O(10 sec)), 
any scenario where the supernova has its energy drained more rapidly is 
ruled out. In this way it is possible || to put a constraint on the majoron 
luminosity. Assuming the p particles are heavy, the largest contribution to 
the majoron luminosity will come from the processes vv — > xXi v ' ~~ * v Xi 
and v — > vx, where the last is matter induced, and is a consequence of the 



difference in energy between helicity states [33]. Choi and SantamariaQ] give 
the constraints from these processes as 



vv 



v' 



; GeV < 2 x 1Q _ 
vi v 1 

v X : > 3 x ID" 5 , (47) 

sec m u 

v -> vx : Z 8 x 10~ 7 . (48) 

The scattering constraint (Eq. |46| ) is slightly stronger than that from nucle- 
osynthesis (Eq. |4lD , whereas the off-diagonal inverse decay constraint (Eq. ^) 
is some two orders of magnitude weaker than the corresponding constraint 
from nucleosynthesis (Eq. |43]), and the diagonal argument (Eq. [48]) has no 
analog for the early Universe. As with the nucleosynthesis constraints, none 
of the supernova constraints are relevant if we confine ourselves to "reason- 
able" values, m u < 30 eV and v > 100 GeV. The constraints were based 
on requiring the majoron luminosity to be less than 3 x 10 53 erg sec -1 for a 
0.85 solar mass core, with a radius of 10 km, at a temperature of 50 MeV. 
Although this condition is somewhat ad-hoc, similar criteria in the study 
of axion emission have proven useful as a guide to more detailed numerical 
studies [ |35]j . 

In the present context of a massive majoron we need to add inverse ma- 
joron decay as an additional production mechanism. Actually, inverse ma- 
joron decay already occurs in standard majoron models. The decay of a 
massless majoron into two positive helicity neutrinos can proceed as the 
analog to the helicity flipping neutrino decay discussed by Choi and Santa- 
maria M. The phase space considerations are similar and the matrix element 
is the same, so we expect a similar constraint. If the gravitationally induced 
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majoron mass is larger than the effective neutrino mass, then the phase space 
will be enhanced. We therefore estimate the constraint from inverse majoron 
decay by scaling to Eq. [IB], 

— * § x 10- 7 -^-, (49) 
V\ 100 kev 

where the neutrino effective mass is of order 100 keV. 

The supernova constraints do not apply if the majorons are so strongly 
coupled that they are trapped inside a "majoron-sphere" which has a black 
body luminosity of < 10 53 erg/ sec. Choi and Santamaria estimate that the 
interaction v\ —> vx wm trap the majorons if 

f^LV > 3 x 10 -3. (50) 
MeV V v x J y ' 

A comparison with the constraints from nucleosynthesis (Eq. [|l]) shows that: 
if majorons are coupled strongly enough to be trapped, then they would 
be brought into equilibrium in the early Universe. Thus, trapping is not 
a viable way of evading the supernova constraint^. Rather, the supernova 
constraints are an extension of the nucleosynthesis constraint. This is a 
general result, the thermal conditions in the supernova are not so different 
than in the early Universe, but the time scale required for trapping is shorter 
than the expansion rate in the early Universe at comparable temperatures. 
For example to trap in the core requires an interaction time less than the light 
crossing time of 3 x 10 -5 sec, whereas the expansion time at T = 50 MeV is 
2 x 10~ 4 sec. Actually, there is an even stronger requirement ]57| : to avoid 



excess heat transport in the supernova core requires significantly larger cross- 
sections than needed just for trapping, which would lead to even stronger 
disagreement with the nucleosynthesis constraint. 

The main weakness, apart from numerical uncertainty, with the con- 
straints summarized above is that they rely on the standard cold core-bounce 
model of Type II supernovae being correct, but that scenario may not be con- 
sistent with the hypothesized majoron parameters. The crucial point is that 
if lepton number is strongly violated than the lepton degeneracy may be 
erased on infall. The resultant increase in entropy may cause a "thermal" 



bounce [38|, and it is unknown whether or not such a scenario can produce 



the observed neutrino pulse from SN 1987A. 



3 



For a possible exception to this, see ref |3£ 
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Other constraints may arise if particles produced in the supernova de- 
cay in flight to observable particles while on the way to Earth. In the ma- 
joron/neutrino sector, the most readily observable particle is the V e . Two 
scenarios have been discussed, a) Heavy neutrinos emitted from the neu- 
trino sphere can decay producing a 10 — 20 MeV V e signal. Such a signal 
would have to be temporally separated from the primary neutrino pulse in 



order to be readily identifiable |p9 |. b) High energy majorons from the core 



can decay producing a ~ 100 MeV signal [40]. Such a signal can be identi 



fied by its spectrum. Indeed, a single high energy event coincident with the 
primary neutrino burst would be significant. Since the neutrino detection 
cross-sections increase with energy, only a fraction of the supernova energy 
must be emitted in high energy particles to produce a signal. The lack of any 
high energy events coincident with SN 1987A therefore places a constraint on 
the production processes for majorons that decay into 77 e 's which is ~ 1 order 



of magnitude more stringent in the amplitude than the limits in Eq. [46-48 



8 Summary 

In this paper we argue that physics at the planck scale may have interesting 
consequences for majoron models. Majoron models contain a global sym- 
metry which, when broken, provides both for neutrino masses and relatively 
rapid neutrino decay; however, gravity need not respect global symmetries. 
Unless the global symmetry arises as an automatic consequence of a local 
symmetry, there is no reason to expect the symmetry to survive at all in the 
low energy phenomenology. To this end we have constructed explicit models 
where B — L is a spontaneously broken local gauge symmetry and the ma- 
joron arises as the Goldstone boson of a spontaneously broken 'accidental' 
U(l) symmetry. By varying the B — L charge of a new scalar field in the 
model, the global symmetry can be protected against operators with dimen- 
sion less than some power d. The most interesting cases are those where 
d > 4, in which case the explicit breaking of lepton number is suppressed by 

The result of explicitly breaking L is that the majoron gets a mass, which 
may have interesting consequences. In addition, by invoking an automatic 
symmetry to protect against explicit symmetry breaking the majoron is nec- 
essarily a linear composition of two scalar fields, only one of which couples to 
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neutrinos. As a result, the majoron coupling to neutrinos may be suppressed 
from its value in simpler majoron models. 

By appropriate choices of parameters the majorons may be stable, and 
thermal relic majorons may play the role of dark matter. Alternatively, they 
may decay into light neutrinos with possible consequences for galaxy forma- 
tion, or be so light and weakly coupled that they have almost no observable 
consequences at all. If the dimension of the gravitational operators is large 
enough, (d > 10), there is the possibility of a coherent energy density in 
majorons, similar to that which arises in axion models. 

If the majoron mass lies in the MeV region there will be consequences 
for big bang nucleosynthesis. For lighter (keV) majorons, as long as neu- 
trino masses lie below 30 eV neither neutrino annihilation into majorons nor 
inverse majoron decay are likely to pose a problem for nucleosynthesis. Su- 
pernova constraints on majoron models are hardly altered by the majoron's 
mass, unless it exceeds a few hundred MeV, which would suppress thermal 
majoron production in supernovae. Suppression of the majoron-neutrino cou- 
pling can soften the constraints from supernovae, even for light majorons. 

Finally, there is not complete freedom in choosing the B—L charges for the 
scalar fields in the theory. An inappropriate choice implies the formation of 
a network of cosmological strings and domain walls that cannot dissipate via 
self-intersection. The energy density in such a network would be prohibitively 
large, and so such choices are proscribed. 
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Figure 1: Leading contribution to g eex . 



Figure 2: Scalar potential showing small explicit breaking. 
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